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ABSTRACT 

C"| , We present panoramic Spitzer/MIPS mid- and far-infrared (MIR/FIR) and GALEX 

Cln' ultraviolet imaging of the the most massive and dynamically active system in the 

local Universe, the Shapley supercluster at z=0.048, covering the five clusters which 
make up the supercluster core. We combine these data with existing spectroscopic 
data from 814 confirmed supercluster members to produce the first study of a local 
rich cluster including both ultraviolet and infrared luminosity functions (LFs). This 
joint analysis allows us to produce a complete census of star-formation (both obscured 
and unobscured), extending down to SFRs~0.02-0.05 M Q yr _1 , and quantify the level 
of obscuration of star formation among cluster galaxies, providing a local benchmark 
for comparison to ongoing and future studies of cluster galaxies at higher redshifts 
' with Spitzer and Herschel. The GALEX near-ultraviolet (NUV) and far-ultraviolet 

, (FUV) luminosity functions (LFs) obtained have steeper faint-end slopes than the 

■^j- ■ local field population, due largely to the contribution of massive, quiescent galaxies at 

MFUV^t— 16. The 24/mi and 70/im galaxy LFs for the Shapley supercluster instead 
have shapes fully consistent with those obtained for the Coma cluster and for the local 
field galaxy population. This apparent lack of environmental dependence for the shape 
of the FIR luminosity function suggests that the bulk of the star-forming galaxies that 
make up the observed cluster infrared LF have been recently accreted from the field 
and have yet to have their star formation activity significantly affected by the cluster 
environment. We estimate a global SFR of 327M Q yr _1 over the whole supercluster 
jJJ ■ core, of which just ~20 per cent is visible directly in the ultraviolet continuum and ^80 

per cent is reprocessed by dust and emitted in the infrared. The level of obscuration 
(Ltr/ Lpuy) in star-forming galaxies is seen to increase linearly with Lk over two 
orders of magnitude in stellar mass. 

Key words: galaxies: active — galaxies: clusters: general — galaxies: evolution — 
galaxies: stellar content — galaxies: clusters: individual (A3558) — galaxies: clusters: 
individual (A3562) — galaxies: clusters: individual (A3556) 
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1 INTRODUCTION 

The conversion of gas into stars is one of the most funda- 
mental astrophysical processes regulating galaxy evolution, 
and hence the star formation rate (SFR) which measures 
the rate at which this conversion occurs, building up the 
stellar mass of the galaxy, is a key observable in galaxy 
evolution studies. The harsh cluster environment is well 
known to affect the star formation activity of the mem- 
ber galaxies, as qua ntified through the star-formation (SF)- 
density relation fe.g. lDressler et all 1985: B alogh et al.ll200ol : 



iLewis et al.1l2002l : iGomez et alj|2003h . While various physi- 
cal mechanisms such as ram-pressure stripping, harassment 
and starvation have been proposed to quench star formation 
in infalling spiral galaxies transforming them into the pas- 
sive lenticulars which dominate cluster cores (for reviews see 
e.g. iBoselli fc Gavazzil [200rj : lHaines et ail 120071 ). the domi- 
nant evolutionary pathway(s) remains unresolved. 

Many of these processes that drive the evolution of 
galaxies also shape the luminosity function (LF), one of the 
most basic and fundamental properties of the galaxy popula- 
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tion l|Benson et al.ll2003h . The LFs of cluster galaxies at var- 
ious wavelengths can hence provide quantitative probes of 
how these dense environments affect the fundamental galaxy 
properties such as the overall mass function, stellar masses 
and SFRs. A key question is whether the LF shows a signifi- 
cant environmantal dependence or is instead universal. This 
remains unclear, with some studies s uggesting relatively lit- 



tle variation with environme n t (e.g. de Propris et al.l 1 19981 : 



IChristlein fc Zabludofil 120031 : iRines fc Gelled I2008I ). while 
other find sigificant differences including brighter char- 
acterisitic luminosities and steeper faint -end slopes (e.g. 
IPopesso et al.l [20051 ; iMercurio et al.ll2006h . which could be 
largely ascribed to the diverse morphological com position of 
cluster and field populations (|de Lapparentll2003t ). However 
the optical luminosities of galaxies depend not only on their 
stellar masses but are also strongly affected by the pres- 
ence of young stellar populations and dust, which can bias 
the optical LFs, making it difficult to reliably interpret the 
data. The near-infrared (NIR) instead is much less sensi- 
tive to the effects of dust or star formation, and hence the 
NIR (e.g. JC-band) LF can be considered a reliable estima- 
tor o f the underlying stellar mass function (Bel l fc de Jond 
l200ll ). The ultraviolet (~2 000A) instead provides a measure 
of the recent SFR over the last 10 8 yr in galaxies, being dom- 
inated by emission from you ng stars of intermediate masses 
(2-5 Mm; iBoselli et al.ll200ll l. and hence the ultraviolet LF 
represents a useful tool to quantify the effects of the cluster 
environment on star formation. 

Ultraviolet radiation from young stars can be strongly 
attenuated by the intervening dust, and indeed the ultra- 
violet emission from spiral galaxies comes predominately 
from those slightly older stars (10 7 — 10 8 yrs) that have mi- 
grated away from their birthplaces in the h ighly obscured 
giant molecular clouds (jCalzetti et al.l [2005h . The amount 
of attenuation depends both on the general properties of 
the galaxy (mass, metallicity, morphology) and the relative 
geometry of the dust and young stars, but for many galax- 
ies more than 90% of the UV photons are absorbed by dust, 
corresponding to A(FUV)>2.5 and hence potentially signf- 
icantly biasing the ultraviolet LF. This energy absorbed by 
the dust is however reprocessed as thermal radiation in the 
mid/far- infrared (8-1000/im), and so observations at these 
wavelengths allow us to quantify this energy a nd thus in- 



fer the amount of obscured star formation (e.g. Kcnnicutt 



1998] ; ICalzetti et~ai1l2007l ; iRieke et al1l2009l ; iKennicutt et al 
2009). Optically-thin infrared and radio continuum emission 



provides an inherently extinction-free estimate of SFRs, and 
so the infrared and radio LFs can be considered complemen- 
tary to the ultraviolet LFs, being sensitive to obscured emis- 
sion missed by the ultraviolet LF, and between them provide 
strong constraints on the global SFRs in cluster galaxies and 
on the impact of cluster related environmental processes on 
star formation. 

The first pioneering ultraviolet observations of cluster 
g alaxies were performed u sing the sounding rocket payloads 
of lSmith fc Cornettl l| 19821 ). The FOCA balloon-borne 40cm- 
diameter telescope surveyed the l ocal Coma and Abell 1367 
clusters l|Donas et al.lll990l . ll99lT ). obtaining the 2000 A ul- 
traviolet luminosity function for Coma. They found that de- 
spite Coma being the archetype of the elliptical-rich and con- 
centrated cluster, in the ultraviolet it is strongly dominated 
by spirals, and its surface density shows low concentration 



and contrast over the field. ICortese et ail (|2003l ) then com- 
bined the FOCA and later FAUST UV cluster surveys to 
produce a composite UV luminosity function for the Virgo, 
Coma and Abell 1367 clusters. This was found to be well 
fitted by a Schechter function with M^ v =— 18.79±0.40 and 
a=— 1.50±0.10 values that did not differ significantly from 
the local UV luminosity function of the field. More recent 
deeper studies of the UV LFs of Coma and Abell 1367 with 
GALEX showed steeper faint-end slopes than for the field, 
due to the co ntribution of passively-evol ving galaxies at faint 
magnitudes jCortese et al.ll2005l . l2008al '). 

The Infrar e d A stronomical Satellite (IRAS; 
Neugebauer et all Il985l ) provided the first opportunity 
to study the infrared properties of cluster galaxies and to 
compare them with the field population. In a pointed IRAS 
survey of the Hercules clust er at z=0.03 6 reach ing 50mJy 
at 60/im (L/ fl ~5xlO 9 I/ ), lYoung et all (| 19841 ) identified 
twenty-four 60/im sources as spiral galaxies belonging to 
the cluster, producing a luminosity function that could be 
fit by a double power law. Notably galaxies classified as 
ellipticals or SOs were absent from the infra red data. ISO 
and m ore recently the Spitzer space telescope (jWerner et al.l 
l2004h have allowed much more sensitive infrared studies of 
star formation in clusters, revealing ubiquitous obscured 
star formation among cluster galaxies whose contribution 
was underestima t ed by >10x by prev ious optical studies 
jDuc et alJ 120021 ; iMetcalfe et~aH 120051 ). For loc al clusters 
the resultant SFRs remain low (^lOMay r" 1 ; iDuc et al.l 
|2002| ; IWolf et all |2009| ; lHaines et all l2010ai ) and the MIR 
luminosity fu nctions are con s istent with those found in 
field regions |Bai et alJ [20061 . 120091 ). At higher redshifts, 
galaxies with obscured SFRs of 30-100 M© yr -1 become 
commonplace, indicative of rapid evolution, but largely 
paral l eling that evolution found in the field (IZheng et al.l 
120071 ; lHaines et all l2009al ; IVulcani et ail |2010| ), consistent 
with the se galaxies represe n ting a recently accreted pop- 
ulation dGeach et al l |2006| . l2009al ; iMarcillac et~aH 120071 : 
iKovama et al.ll2010l). 

In lMerluzzi et alJ (|2010l . Paper I) we introduced the AC- 
CESS multi-wavelength survey of the Shapley supercluster 
core (SSC) as well as presented an analysis of the K-band 
data. The target has been chosen since the most dramatic 
effects of environment on galaxy evolution should occur in 
superclusters, where the infall and encounter velocities of 
galaxies are greatest (>1 OOOkms" 1 ), groups and clusters 
are still merging, and significant numbers of galaxies will be 
encountering the dense intra-cluster medium (ICM) of the 
supercluster environment for the first time. Moreover, the in- 
terplay of a variety of physical processes need to be invoked 
to d escribe the ov erdensity of baryons on supercluster scales 
(e.g. lFabianlll99ll ), which emphasises the cosmological con- 
text of studies involving the effects of star formation, AGN 
and associated gas processes. 

The Shapley supercluster was identified as the largest 
overdensity of galaxies and clusters in the z<0.1 Universe in 
an APM survey of ga laxies covering the entire southern sky 
dRavc haudhurv 1989]) and confirmed by ear ly X-ray studies 
l|Ravchaudhurv et al.lll99ll;|Pav et al.lll99ll). In the original 
Abell catalogue (|Abell. Corwin fc Olowin 19891 ). the core of 
the Shapley supercluster was found to be represented by a 
single cluster (A3558=Shapley 8) of optical richness 4, in- 
consistent with the observed velocity dispersion of the clus- 
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Figure 1. The multi- wavelength photometric coverage of the ACCESS survey. The shaded region indicates the WFI B, R - band im aging 
of the Shapley Optical Survey. Coloured contours represent the surface density of R<21 galaxies (Fig. 1 of lHaines et al.l l l2006af )). The 
respective coverages of the GALEX FUV/NUV (blue circles), WFC AM/ftf-band (red dashed-line), Spitzer 24^tm (magenta lines) and 
70/im (green lines) imaging are all indicated. 



ter. Early X-ray observations from Einstein and ROSAT 
l|Breen et al.lll994l ) resolved this by discovering that the core 
in fact comprises of a chain three merging clusters, A3558, 
SC1327-312 and SC1329-313. The SSC in this paper refers 
to these three clusters, along with the two associated Abell 
clust ers on either side, A3562 and A3556 (see lMercurio et"al] 
2006, for a general description). 

In this paper we attempt to build a complete and un- 
biased census of star formation in galaxies belonging to the 
SSC using a multi-wavelength analysis combining GALEX 
ultraviolet and Spitzer infrared observations, in order to 
circumvent any biases that may exist in previous surveys 
in which star formation rates in cluster galaxies were esti- 
mated from a single SFR indicator. By assembling a com- 
prehensive pan-chromatic (FUV-FIR) dataset for the SSC 
we are able to fully account for both obscured and unob- 
scured star-formation, and by dealing with global photomet- 
ric measurements, we are independ ent of the aperture bi- 
ases (IKewlev. Jansen fe Gellerll2005l ) prevalent in many pre- 
vious studies based on optical spectra obtained from fibre- 
fed spectrographs (e.g. iLewis et~al 1 120021 ; iGomez et al1l2003l ; 
lHaines et al.ll2007h . In §[2] we present the ultraviolet and new 
infrared datasets and describe their reduction and resultant 
source catalogues. In §[3] and §[4] we present the ultraviolet 
(NUV/FUV) and infrared (24pm, 70/xm) galaxy luminosity 
functions for the Shapley supercluster. We then perform a 
joint analysis of the ultraviolet and infrared data to measure 
their relative contributions to the global supercluster SFR 



budget in § [S] and quantify the level of obsuration among 
cluster galaxies (both individually and globally), before pre- 
senting the discussion in §[6] 

In a companion paper l|Haines et al.ll2010bl . Paper III) 
we examine in detail the nature of star formation in galaxies 
within the Shapley supercluster, finding a robust bimodality 
in the /24//K colours that can be understood as the well- 
known split into star-forming and passive galaxies, as well 
as taking adva ntage o f avail able 1.4 GHz VLA radio contin- 
uum data from lMillerl (|2005f ) to show via examination of the 
panchromatic (FUV-FIR) SEDs that the bulk of the global 
supercluster SFR comes from quiescent star formation in 
normal infalling spiral disks who have yet to be affected by 
the supercluster environment. In a further follow-up paper 
(Haines et al. 2010c, in preparation) we will examine the 
environmental trends in star formation, and identify and 
quantify various classes of galaxies in the process of being 
transformed from star-forming spiral to passive SO. 

This work is carried out in th e framework of the 
joint research programme ACCES^f] (|Merluzzi et al. I l2010l . 
|http://www.na.astro. it/ A CCESS ) aimed at distinguishing 
among the mechanisms which drive galaxy evolution across 



4 A Complete Census of Star-formation and nuclear activity in 
the Shapley supercluster, PI: P. Merluzzi, a collaboration among 
the Universities of Durham and Birmingham (UK), the Italian 
National Institute of Astrophysics' Osservatorio di Capodimonte 
and the Australian National University (ANU). 



4 Haines et al. 



different ranges of mass by their interactions with the en- 
vironment. When necessary we assume Qm=0.3, Qa~0.7 
and Ho=70 kms _1 Mpc _1 , such that at the distance of the 
Shapley supercluster 1 arcsec is equivalent to 0.96 kpc. All 
magnitudes are quoted in the Vega system unless otherwise 
stated. 



2 DATA 

The ultraviolet and mid-infrared data analysed in this pa- 
per are complemented by existing panoramic optical B- 
and fl-band imaging from the Shapley Opti cal Survey 
(SOS: iMercurio et al.ll2006l; lHaines et al l [20 06a ) and near- 
infrared A-band imaging i|Merluzzi et al.ll2010l ), which cover 
the same region. Figure \T\ shows the Shapley superclus- 
ter core (indicated by isodensity contours) with superposed 
the relative coverages of the different aspects of our multi- 
wavelength survey. The ultraviolet (blue circles) and mid- 
infrared (24/im magenta; 70pim green) photometry cover es- 
sentially the same regions, the bulk of which are also covered 
by the A-band imaging (red dashed-lines) , while the ear- 
lier SOS covers a slightly narrower region. We thus have full 
multi- wavelength coverage of the main filamentary structure 
connecting A3562 and A3558, and the extension to A3556. 



2.1 The Shapley Optical Survey 

The SOS comprises wide-field optical imaging of a 2.2 deg 2 
region covering the whole of the Shapley supercluster core. 
The observations were made from March 2002 to April 2003 
using the Wide Field Imager (WFI) camera, on the 2.2m 
MPG/ESO telescope at La Silla. The SOS is made up of 
eight contiguous WFI fields of 34 x 33 arcmin 2 (pixel scales 
of 0.238 arcsec), each with total exposure times of 1500s in 
B and 1200s in R, and typical FWHMs of 0.7-1.0 arcsec, re- 
sulting in galaxy catalogues that are both complete and reli- 
able to fl=22.0 and £?=22.5. Full details of the observations, 
data reduction, and the production of the galaxy catalogue 
are described in IMercurio et ah! l|2006l ). The high quality of 
the optical images has allowed morphological classifications 
and structural para meters to be derived fo r all supercluster 
galaxies to ~A'/^+3 dGargiulo et al.ll2009l) within the auto- 
mated 2dphot environment ( La Barbera et all [2008). 



2.2 A-band and spectroscopic data 

The A-band survey of the SSC was carried out with the 
Wide Field infrared CAMera (WFCAM) on the 3.8m United 
Kingdom Infra-Red Telescope (UKIRT) in April 2007 as 
part of the UKIRT service programme (P.I. R.J. Smith). 
The WFCAM data consist of five complete 52x52arcmin 2 
tiles (made up of four interleaved multiframe exposures), 
covering a region of 3.0 deg 2 , with exposure times of 300s, 
0.4arcec pixel scale, FWHMs of 0. 9-1.2 arcsec, and rea ch 
A=19.5 at 5a. For more details see lMerluzzi et all ([2010). 

The Shapley supercluster core has been observed 
by a variety of redsh ift su rv eys (lOuintana et alj 
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1995I; iBardelli et all Il998l. |2000|; iKaldare et alj 2003 



Drinkwater et alj|2004 ISmith et al.ll2007l ; ICava et alj|200£ ; 
Jones et al.ll2009l ) resulting in 964 galaxies in our A-band 
catalogue having redshift information, of which 814 lie in 



Figure 2. The completeness of the GALEX NUV (thick red) 
and FUV (thick blue) photometric catalogues as a function of 
rriAB-, averaged over the four GALEX images. The thin black 
lines indicate the completeness levels for the individual GALEX 
NUV and FUV images. 



the velocity range of the SSC (10 700<cz<17000 kms" 1 ). 
Of these, 415 are from the AAOmega-based spectroscopic 
survey of ISmith et"aT] l|2007f ). We have redshifts for all 107 
A<12.3 (A<A*+0.6) galaxies within the A-band survey. 
The 90% and 50% spectroscopic completeness limits are 
A=13.25 (AT* +1.5) and A=14.3 (A* +2.5) respectively. 



2.3 The GALEX observations 

The UV data analysed in this paper consist of four con- 
tiguous GALEX pointings (blue circles in Fig. [1} making 
up the Cycle 4 Guest Investigation proposal GI4-098 (P.I. 
R.J. Smith), and which were observed in May 2008. The 
GALEX field of view is circular with diameter 1.2°, with 
pixel scale 1.5 arcsec. Each pointing is simulaneously ob- 
served in both the far-UV (FUV) and near-UV (NUV) bands 
with effective wavelengths of 1516A and 2267A, and spatial 
resolutions of 4.3 and 5.3 arcsec respectively. The exposure 
times are in the range 1238-1688 s, i. e. comparable to the 
GA LEX Medium Imaging Survey . See iMartin et al.1 (|2005l ) 
and iMorrissev et alj l|2005l . [2007 ) for details regarding the 
GALEX instruments and mission. 

Sources were detected and measured from the GALEX 
images using SExtractor (|Bertin fc Arnoutsll996t ). As the 
NUV images are deeper than the FUV images (in terms 
of source counts), we used the NUV images for detection 
and measured the FUV flux in the same ape rture as for 
the N UV. Following IMercurio et al.1 (|2006l ) and iGrav et al.1 
(2009), we adopted an optimized, dual ('hot' and 'cold') con- 
figuration, in order to deblend high-surface brightness ob- 
jects that are close on the sky in projection, and simultane- 
ously avoid the spurious shredding of highly structured spi- 
ral galaxies into multiple star-forming regions. The FUV and 
NUV magnitudes are corrected for Galacti c extinction us- 
ing the lSchlegel. Finkbeiner fc Davis] (|l998i ) reddening map 
an d the parametrization of the Galact ic extinction law given 
bv lCardelli. Clayton fc Mathisl l| 1989 1 ) for a total-to-selective 
extinction ratio of fly =3.1, resulting in the conversion fac- 
tors A FUV =7-9E(B - V) and A N uv=8.0E(B - V). For 
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Figure 3. Completeness function of our 24/im catalogue. Each 
one of the five black histograms corresponds to the mean of the 
two simulations for each field. The red histogram is the average 
over the whole survey area. According to these simulations, the 
survey is 90% complete to 0.35mjy. 



our fields E(B — V) — 0.048-0.058, resulting in corrections 
of ~0.40 mag. 

The flux detection limits and completeness of each 
GALEX image were determined by individually insert- 
ing 500 simulated sources at random positions across the 
GALEX image for a range of fluxes within a given magnitude 
bin, and determining their detection rate and recovered mag- 
nitudes, using identical extraction procedures. This is then 
repeated for all magnitude bins of interest. As for all but 
the brightest galaxies, sources are unresolved in the GALEX 
images, we used the NUV and FUV point-spread functions 
provided by the GALEX website as our simulated sources. 
The resulting completeness levels as a function of tuab for 
the NUV (red curve) and FUV (blue curve) over the range 
19.5<m,As<23.5 are shown in Fig. [2] We expect our UV 
catalogues to be 90% complete to rriAB{NUV)=22.0 and 
mAB{FUV)=22.5. Based on the ultraviolet c omponent of 
the co mbined FUV+24/xm calibration (Eq.[T]) of lLerov et al.l 
(l2008t). whi c h itse lf comes from the FUV-SFR relation of 
ISalim et al.l l|2007h the FUV completeness limit corresponds 
to a SFR of 0.014 M yr _1 at the SSC distance. 

The UV detections were cross-matched with the optical- 
NIR catalogues using a search radius of 5 arcsec (except 
for two large edge-on spirals for which a search radius of 
15 arcsec was required). For the remainder of this article, we 
limit our analysis of the GALEX data to either the 2.83 deg 2 
region also covered by if-band photometry, or the 2.26 deg 2 
region also covered by the SOS B- and i?-band photometry. 



2.4 The Spitzer observations 

The panoramic Spitzer mid-infrared observations of the 
Shapley supercluster core were carried out over 27-30 Au- 
gust 2008 within the Cycle 5 GO programme 50510 (PI: 
OP. Haines). The observatio ns consist of five contiguous 
mosaics observed with MIPS i Rieko el al. 21111 ll ) in medium 
scan mode, with coverages at 24/im and 70^m shown as ma- 
genta and green boxes respectively in Figure \T\ This scan 
rate leads to 10 exposures each of 4 sec, for each pixel along 
a single scan. As the usable 70/im array only covers half the 
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Figure 4. Comparison of the Spitzer survey depths at 24/xm (blue 
bar) and 70/^m (green bar) to the luminosity-dependent infrared 
SEDs of iRieke et all (2009) redshifted to the Shapley superclus- 
ter. The response functions of the Spitzer/MYPS filters are indi- 
cated by the filled curves. 



full MIPS array width, we used a scan leg spacing equal to 
the half-array width (160 arcsec). This strategy results in a 
homogeneous coverage across both 24/im and 70/im mosaics, 
with each point covered by two 24/im scans and one 70/im 
scan during a single observation. 

The reductio n of the MIPS data was perf ormed with the 
MOPEX package (|Makovoz fc Mar lcau 2005, version 18.1.5) 
and with th e Germanium Repr ocessing Tools (GeRT, ver- 
sion 060415; iGordon et aHl200 5f|. 



2.4.1 24/im data 

We reduced each of the five fields (corresponding to the five 
AORs) separately, with the aim to allow for cross-checking 
during the reduction. Since the mosaics delivered by the 
pipeline were affected by dark latents along the scan direc- 
tion, we started by applying to the Basic Calibrated Data 
(BCD) the so-called 'self-calibration', consisting of dividing 
the BCDs by a flat-field derived from the normalised median 
of all BCDs in each AOR (using the MOPEX module flat- 
field. pi). We then produced a mosaic for each field usi ng the 
MOPEX module mosaic.pl l|Makovoz fc Khanl I2OO5T ) . The 
average effective exposure time per pixel in the final mo- 
saics is 84s. 

For the detection of the sources, MOPEX relies on 
an empirical determination of the Point Response Function 
(PRF) from the actual data. To derive the PRFs for our 
mosaics, we first produced a theoretical PRF using STiny- 
Tim, a program for computing PRF models starting from 



This software, as well as full information on data anal- 
ysis is available at the Spitzer Science Center website 

1 http:/ / ssc. spitzer. caltecKed u/ dataanalysi stools'/\ 
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Figure 5. Relative error of the measured 24/j.m flux as a function 
of the input flux, as derived from simulations. Points correspond 
to the simulated point sources. The thicker red curve corresponds 
to the mean error, and the top and bottom curves to ±lcr from the 
mean error. These curves are used to estimate our flux measure- 
ment errors. Mean errors increase for the lower fluxes because 
faint sources are preferentially detected on the top of po sitive 
noise fluctuations, as already noticed bv lFadda et al.l (|2006h ■ 
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Figure 6. The relation between the x 2 goodness-of-fit value of 
the MIPS point response function to the 24/mi emission and the 
total 24/im flux for Shapley supercluster members (open sym- 
bols) and stars (red stars). The blue dotted curve indicates the 
lower limit for the x 2 distribution of 24^tm supercluster galax- 
ies clearly resolved by MIPS. Those supercluster galaxies whose 
24/^m emission is unresolved are indicated by shaded symbols. 



instrumental parameters. The model PRF was then used to 
extract a first catalogue of point sources. From the 30-40 
brightest sources of each field, we selected ~10 'bona- fide' 
point sources. These were then used to derive the PRF with 
the MOPEX module prf_estimate.pl. The catalogues were 
extracted using the MOPEX module apex_lframe.pl, which, 
besides source detection, also derives fluxes by fitting the 
PRF to the detected sources. The fit is performed inside a 
normalisation radius, which we fixed to 5.9 arcsec, which in- 
cludes the central, brightest part of the PRF. In total 25 718 
sources were detected with SNR>3, which was the minimum 
signal -to-noise ratio w e set for a potentially real detection 
(c.f. iFadda et alj|2006l fl 

To assess the reliability of source extraction and flux 
estimates and to determine the completeness of our cata- 
logues, we repeated exactly the same process as outlined 
above on a set of simulated images. To create the simulated 
images, we first subtracted from the original mosaics the cat- 
alogued sources, using the MOPEX module apex_qa.pl. On 
these 'empty' fields, we then added the same point sources 
of the original catalogues, in random positions, but avoiding 
the areas previously occupied by the 'real' sources, where 
residuals of the subtraction could be present. To create the 
new point sources, we multiplied our PRFs by the fluxes in 
the catalogues. This process was conceived in order to put 
the artificial (but realistic) point sources on as realistic a 
background as possible. We repeated the process twice on 
each of the five fields. The simulations were performed using 
FORTRAN codes plus the apex_qa.pl module. 

We computed the completeness function by dividing, 
for each flux bin, the number of detected sources by the 
number of input sources. The completeness functions of our 



6 Notice that our SNR values are estimated with a different al- 
gorith m of MOPEX, not yet available at the time of lFadda et al.l 
(2006) . The new S NR v alues should be doubled to be compared 
with lFadda et al.l ||2006TI 



24/xm data are given in Fig. [3] where we show the results for 
the five fields and their average, from which it is apparent 
that the completeness functions are fully consistent among 
each other. From this figure we conclude that our 24^im data 
are 90% complete to a flux of 0.35 mjy, which we adopt as 
our survey limit. In Figure [4] we show how this limit relate 
to the luminosity-dependent infrared SEDs of iRieke et al.l 
(2009) placed at the distance of the Shapley superclus- 
ter. From this figure we see the 24/xm limit corresponds to 
logL/ fl (ergs~ 1 )=42.46 or L IR =7.5x 10 8 L Q . Based on the 
in frared component o f the combined FUV+24/im calibration 



Lerov et al 



Ricke et al 



2008) , or equivalently the 24/^m calibration 

2009) , this corresponds to an obscured SFR 
i.e. comparable, but marginally higher than 



of 0.05M Q yr" 
the FUV-based limit. 

Fig. [5] shows the relative error of the measured flux as 
a function of the 'true' (input) flux as derived from simula- 
tions. Each point corresponds to a point source, while the 
red curves correspond to the mean error (the central, thicker 
curve) and to the 1-a limits (top and bottom curves). These 
last curves were used to estimate our flux measurement er- 
rors. These errors could then added to the uncertainty of 
~4% on the absolute flux calibration at 24/im (see the MIPS 
Instrument Handbook at the Spitzer Science Center web 
site) and to the 2% uncertainty in the aperture corrections 
(see below) to obtain the total uncertainty on the fluxes. 

Since the fluxes were estimated inside a radius 5.9 arc- 
sec, we must apply aperture corrections to obtain the total 
fluxes from the sources. To this purpose, we derived growth 
curves for a number of different cases: (i) model point sources 
with black-body temperatures of 15K, 30K, 50K, 500K (cre- 
ated with STinyTim); (ii) the growth curve given in the 
MIPS Data Handbook (Fig. 3.2 in versio n 3.3.1); (iii) the 
growth curves given bv lFadda et"afl (120061 ); (iv) the growth 
curve in the Spitzer Science Center web site. We also mea- 
sured the growth curve from one bright point source at the 
centre of a field, which gave a result consistent with the oth- 
ers, but with a small systematic positive offset, which we as- 
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Figure 7. Completeness function of our 70/tm catalogue. The 
six black histograms correspond to the results of the six simula- 
tions, while the red histogram is their average. According to these 
simulations, the survey is 90% complete to 25mjy. 



cribe to incomplete subtraction of the surrounding sources. 
All of these growth curves agree, to within 2%, to determine 
an aperture correction of 1.69±0.03 for our 24/im fluxes. 

In total we associate 4663 /24>0.35mJy sources with 
galaxies in the i\-band image, of which 490 have confirmed 
redshifts, 394 of which place them in the Shapley superclus- 
ter. We have redshifts for all /24>26mJy sources, and 50% 
redshift completeness at ^5 mjy. 



2.4.2 2Jt\im photometry of extended sources 

The vast majority of sources detected in our 24/im mo- 
saics are unresolved at the spatial resolution of the MIPS 
instrument (FWHM~ 5"), and so the optimal estimate of 
their 24/im fluxes should be obtained via the point-spread 
function fitting approach of MOPEX. However, at the rela- 
tively low redshift of the Shapley supercluster, the brightest 
galaxies have radii somewhat larger (5-20") than the MIPS 
point-spread function, and can be seen to be resolved by the 
24/im imaging. In these cases, the sources are poorly fitted 
by point-source functions in the 24/im images, as witnessed 
by the large residuals left after subtraction of the best-fit 
point source from the images, and the resultant fluxes pro- 
duced by MOPEX will be systematically underestimated. 

We can separate those galaxies resolved by MIPS, and 
those for which the galaxies appear as point-sources, using 
the x 2 goodness-of-fit value produced by MOPEX when fit- 
ting the 24/tm emission by the MIPS point spread function. 
Figure [6] shows the relation between this \ 2 goodness-of-fit 
value and the total 24/mi flux for supercluster galaxies (open 
symbols) and stars (red stars). For /24>5mJy we can dis- 
tinguish two sequences, one corresponding to galaxies and 
a second with lower \ 2 values due to point sources of stel- 
lar origin. Along this second sequence, we also find a small 
population of supercluster galaxies (shaded circles) whose 
X 2 values are closer to those of stars having the same 24/im 
flux levels, indicating that their 24/im emission is unresolved 
by MIPS. 

To deal wit h these extended sourc es, we also ap- 
plied SExtractor (|Bertin k, Arnoutsl Il996t ) to the median- 
subtracted 24/im mosaics (fixing the sky background level 




0) 



-0.5 



20 



50 



100 



200 



Flux TnnJyl 



Figure 8. Relative error of the measured 70/tm flux as a function 
of the input flux, as derived from simulations. Points corresponds 
to the simulated point sources. The thicker red curve corresponds 
to the mean error, and the top and bottom curves to ± la from 
the mean error. These curve are used to estimate our flux mea- 
surement errors. 



in SExtractor to zero), performing photometry in a range 
of seven circular apertures with diameters in the range 
10-90 arcsec, and correcting the apertures fluxes based on 
growth curves obtained in § 12.4.11 For bright galaxies (i?<15 
or A"<12.5) and those galaxies for which the 24/im pho- 
tometry was poorly fit by a point-spread function within 
MOPEX, we then compared the flux estimate within the 
30 arcsec diameter aperture obtained by SExtractor with the 
MOPEX estimate, prefering the SExtractor-based flux esti- 
mate in cases where it was found to be >3<r higher than 
that from MOPEX. If this were the case, we then consid- 
ered the larger apertures in turn, again preferring the flux 
estimate from the larger aperture if it were found to be >3er 
higher than that of the smaller aperture. In total, of the 
25 718 sources detected by MOPEX in the 24/im mosaics, 
we identified 231 sources for which the aperture-corrected 
SExTRACTOR-based flux estimate was used instead of the 
MOPEX one. However, the fraction of galaxies for which 
the 24/im flux estimate is SExTRACTOR-based increases with 
24/im flux, such that they represent the majority of galaxies 
with /24>2.5mJy. 



2.4.3 70fj,m data 

The brightest sources in the 70/im images were affected by 
negative side-lobes due to the online filtering when process- 
ing of the BCDs with the automated pipeline. To solve this 
problem, we performed a double-pass filtering using the Ger- 
manium Reprocessing Tools. The double-pass filtering con- 
sists of an initial filtering along the columns of the BCDs, 
followed by a high-pass time median filtering, avoiding the 
pixels near to bright sources, which for this purpose are pre- 
viously masked. This process minimizes data artifacts (al- 
lowing us t o eliminate the sid e-lobes) while preserving the 
calibration jFraver et al.ll2006T l. 

A single mosaic was created starting from all of the 
~9700 filtered BCDs. The subsequent processing was anal- 
ogous to the 24/im data. The average exposure time per 
pixel on the final mosaic turned out to be 42s. The nor- 
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RA 


Dec 


cz 


K 


R 


B-R 


hi 


fro 


FUV 


NUV 


(J2000) 


(J2000) 


(km/s) 


(mag) 


(mag) 




(mjy) 


(mJy) 


(AB mag) 


(AB mag) 


13:23:28.949 


-31:43:12.96 


14036 




18.845 


0.447 


0.313±0.107 


O.OiO.O 


20.290 


20.355 


13:23:43.990 


-31:47:00.55 


14325 




16.540 


0.875 


1.518±0.111 


O.OiO.O 


19.084 


19.726 


13:23:49.351 


-31:42:56.89 


12922 




16.698 


0.987 


2.902±0.112 


19.8±6.4 


21.808 


24.114 


13:24:05.597 


-31:36:36.11 


14711 


16.561 


18.952 


1.145 


O.OOOiO.OOO 


O.OiO.O 


99.000 


99.000 


13:24:20.441 


-31:23:31.05 


13434 


15.203 






2.478±0.112 


34.0±6.6 


20.160 


20.966 



Table 1. Spcctroscopically confirmed supercluster members covered by either our i<"-band or optical imaging. Non-detections in a given 
passband are indicated by a magnitude of 99 or a flux of zero. No value is given in the case when a galaxy is not covered by a particular 
passband. The full version of this table is given in the electronic version of the journal - see Supporting information. A portion is shown 
here for guidance regarding its form and content. 



malisation radius for the PRF fitting was in this case 17 
arcsec, corresponding to the beginning of the first Airy 
minimum for the 70[im PRF. A total of 1740 sources 
were selects with SNR>3. Simulated images were cre- 
ated in the same way as for the 24/im data, except that 
in this case we created six images of the same field. 
The resulting completeness function is given in Fig. 
which shows that we reach the 90% completeness limit at 
25 mjy, which we a dopt as our surve y limit. Based on 
the infrared SEDs of iRieke et all (|2009l ). this corresponds 
to log(L Jfl /L )=9.76 or log(L Jfl fergs- 1 1)= 43.34 (Fig. B 
Based on the direct 70/im SFR calibration of lCalzetti et al.l 
l|2010l ). SFR 7 o M m=I/7o(ergs- 1 )/1.7xl0 43 , this corresponds 
to a SFR of 0.4Moyr -1 , a factor 10 less sensitive than the 
24/xm limit . Alte rnatively, the infrared SFR calibration of 
iBuat et all (120081 ). SFRj«M s yT _1 =(l-7j)10 _9 - 9T (Lj«/£ s ), 
produces an identical SFR limit in th e case when r\ — 0.3 
(following Iglesias-Para mo~et al.l [20041 ). where rj is the frac- 
tion of dust emission due to heating from evolved stars. 

Figure [8] shows the relative flux errors as a function of 
flux, along with the curves encompassing la deviation from 
the mean error, which at our survey limit is 6.5 mjy. To 
this error we added another 7% due to the absolute flux 
calibration at 70/im (MIPS Instrument Handbook) and the 
3% uncertainty in the aperture corrections (see below) to 
obtain the total uncertainty on the fluxes. 

As for the 24pm data, to derive the aperture corrections 
we analysed different growth curves, in particular: (i) four 
model point-sources, created with STinyTim, with black- 
body temperatures of 15K, 30K, 50K and 500K; (ii) three 
growth curves given by the Spitzer Science Center web-site 
(two black-bodies at 15K and 3000K and one power-law 
spectrum). From these curves we derived an aperture cor- 
rection of 1.74, consistent within 3% among all curves. 

We identify 728 sources with /7o>25mJy, of which 589 
are covered by our 7\-band imaging. Of these, 173 have no 
_K"-band counterpart within 5 arcsec, but we expect most of 
these to be either asteroids or high-redshfft ULIRGs/QSOs, 
the former lacking if-band counterparts due to their move- 
ment across the sky in the time between the two sets of 
observations. A further 35 sources were classified as stars in 
the _K"-band image, but we note that almost all of these are 
probably high-redshift AGN rather than stars given their 
high /7o//if ratios. In total, we associate 381 /7o>25mJy 
sources with galaxies in the 7\-band image, of which 194 
have confirmed redshifts, 160 of which place them in the 
Shapley supercluster. 

In Table [1] we present the positions, redshifts, UV, op- 



tical and NIR magnitudes, and Spitzer/MTPS fluxes of all 
galaxies spectroscopically confirmed as belonging to the SSC 
covered by either our K-band or optical imaging. 



3 UV LUMINOSITY FUNCTIONS 

To determine the UV luminosity function of galaxies in 
a cluster a reliable measure of the contribution from 
back/foreground galaxies to the UV counts is required. 
Where spectroscopic information is lacking, we must rely on 
statistical methods to estimate the background contamina- 
tion per magnitude bin, and/or use the UV-optical colours 
to separate supercluster and background galaxies. 

To understand where and how this applies to our anal- 
yses, we show in Figure [9] the UV-optical colour- magnitude 
diagrams of galaxies in the Shapley supercluster core, in 
which the large green (small blue) symbols indicate galaxies 
with (without) redshift information placing them within the 
Shapley supercluster (10 700<cz<17 000 kms -1 ). In each 
panel we see the well known colour bimodality into the 
quiescent red sequence (NUV— i?~6; F UV—R~7) and blu e 
cloud (l<NUV-i?<4) populations (c.f. lHaines et al.ll2008l ). 
As might be expected, the brightest supercluster galaxies in 
the FUV and NUV all have colours indicative of actively 
star-forming galaxies. 

In total 641 galaxies with both UV and 7\-band pho- 
tometry have redshifts, of which 552 are supercluster mem- 
bers. At bright magnitudes (rriAB<18), the estimation of 
the UV luminosity function of galaxies in the SSC can 
largely be performed spectroscopically, since to this point 
our redshift completeness is ^90%, falling to ~50% over 
18<mAs<20. At fainter magnitudes, it is notable that while 
we have redshift information for the vast majority of the 
quiescent galaxy population {NUV— _R>4.5) down to (and 
beyond) our UV completeness limits, for the star-forming 
"blue cloud" populations our redshift completeness drops 
rapidly beyond A/W~19.5 and fW~20. This difference is 
due to the redshift surveys being optically-selected (Bj or R 
bands) and the substantial range in UV-optical colours be- 
tween star-forming and quiescent populations. It will hence 
be necessary to quantify how many of the blue galaxies at 
faint UV magnitudes belong to the Shapley supercluster, 
and how many are background galaxies. Given that in each 
UV magnitude bin, the dominant contribution to the num- 
ber counts will come from this blue star-forming population, 
this is the critical issue in determining the faint-end slope. 
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Figure 9. UV-optical colour-magntiude diagrams of galaxies in the Shapley supercluster core. Large green symbols indicate those 
galaxies spectroscopically confirmed as belonging to the supercluster (10700 < vel < 17000 kms -1 ), blue symbols indicate those galaxies 
photometrically selected as supercluster members. Black dots indicate sources photometrically identified as background galaxies. 



Following ICortese et all l|2008al h we consider two different 
methods to estimate the UV galaxy luminosity function. 



3.1 LF estimation via statistical subtraction of 
field galaxies 

In the first method we statistically account for field galax- 
ies, per bin of UV magnitude, that are expected to lie 
within the area of interest by refe rring to the UV galaxy 
counts published bv lXu et all (120051 ). By matching GALEX 
photometry from 35 MIS (Medium Imaging Survey) and 
3 DIS (Deep Imaging Survey) pointings with optical star- 
galaxy classifications from the SDSS, they produced num- 
ber counts of UV galaxies as a function of NUV and FUV 
magnitude over 24deg 2 to rriAB~22.6 and 1.7 deg 2 over 
22.6<?tiab<23.6. We thus statistically account for contam- 
ination of background galaxies as a function of UV magni- 
tude in bins of width 0.4 mag over 14<rriAs<20 and 0.2 mag 
over 20<rriAs<23.6 (matching the magnitude bins provided 
by Xu et al. 2005), and for each galaxy in a given bin esti- 
mate the probability that it belongs to the Shapley super- 
cluster as pi — (Nsc — Nf)/Nsc, where Nsc is the number 
of UV galaxies in the SSC survey region in that magnitude 
bin, and Nf is t he corresponding galaxy number count from 
IXu et ail (|2005h normalized to the Shapley UV survey area. 
Summing all of the pi within a given magnitude bin, we re- 
obtain the statistical estimate for the number of supercluster 
galaxies in that bin, such that pi can be considered to be the 
contribution of each galaxy to the UV LF. For those galaxies 
with available redshifts, Pi = 1 for 10 700<cz<17000kms _1 
or otherwise. 

We determine the UV k-corrections for each galaxy on 
the basis of its FUV, NUV and (when available) B pho- 
tometry. For star-forming galaxies with FUV— NUV^0.5 
(or FxcxA -1 ), these k-corrections are negligible (<0.05 mag), 
but for quiescent galaxies with FUV— iV[/V~l-3, these can 



reach 0.2-0.3 mag at 2=0.048 (see Fig. 2 of iTrever et ail 
l2005h . Note that rather than fit model spectral energy dis- 
tributions (SEDs) to the galaxy colours to estimate the UV 
k-correction, we simply fit power laws through the two ad- 
jacent bandpasses (FaocA' 3 ), which is generally a reasonable 
description of a star-forming galaxy's SED in the UV. We 
take the local spectral index for the FUV band (from which 
we derive the k-correction) from the spectral index corre- 
sponding to the FUV-NUV colour, and in the case of NUV 
take the local spectral index to be the average of that ob- 
tained from the FUV-NUV and NUV-B colours (if the latter 
is available). 

The primary source of uncertainty (and one which is 
often neglected) in estimates of cluster galaxy LFs is often 
the effect of clustering and cosmic variance on the back- 
ground faint galaxy counts. This effect should be somewhat 
lower in the ultraviolet than in optical or infrared surveys of 
comparable depths, as ultraviolet-selected galaxies are less 
clustered than samples selected at longer wave l ength s, and 
also tend to lie at lower redshifts. iMilliard et al.l (|2007l ) show 
that the median redshift of NUV <22 galaxies is 0.25, with 
little contribution from galaxies at z>0.4, hence any large- 
scale structures at high redshifts should have little or no 
effect on our ultraviolet galaxy counts. The contribution of 
cosmic variance to the uncertainty in galaxy counts within 
a survey region can be estimated from the angular cor- 
relation function w(0) as a 2 = nfl + n 2 J dflidfl2w(6i2), 
where n is the number density of galaxies, and dfii and dD.2 
are elements of the solid angle Q. separated by an angle 6\2 
|Ellis fc Bland-Hawthornll2007l 'l. We include the effects of 
large-scale structure in the uncertainty in background galaxy 
counts in each magnitude bin us ing the the NUV and FUV 
angular correlation functions of IMilliard et al. (2007). 

The FUV (blue symbols) and NUV (red symbols) 
galaxy luminosity functions for the whole 2.8334 deg 2 re- 
gion covered by both GALEX and WFCAM K-bsmd 
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Figure 10. UV-optical colour selection criteria to identify supercluster members, (left panel) NUV — R/B — R UV-optical colours, (right 
panel) B—R/R colour magnitude diagram. Large green symbols indicate galaxies spectroscopically confirmed to belong to the supercluster 
(10700 < vel < 17000 kms -1 ), while red symbols correspond to galaxies known to lie behind the supercluster (vel > 17000 kms -1 ). 
Blue dots indicate galaxies (R < 20 in the left panel) without redshifts. The red lines indicate the selection criteria used to identify 
supercluster members based on (i) their UV-optical colours (left panel) and (ii) their position with respect to the cluster red sequence 
(right panel) 



imaging are presented as open symbols in Fig. 1111 
extending to Mnuv=— 14 and Myw = — 13.5. We fit 
single Schechter functions (shown as dashed curves) 
via a x 2 - m inimization routine, obtaining best-fit values 
M Fuv=-18-94±°;^ and ajw=-1.67±0.03 for the FUV 
LF, and Mntjv=-19.49± u ;3i and a N uv=-1.67±0.035 for 
the NUV LF. The single Schechter fit provides a good 
description for both the FUV and NUV LFs down to 
M NUV =-14 and M F uv = -13.5, with x£~0.9 for 10 dof. 

3.2 LF estimation via UV-optical colour selection 

In the second method we use UV-optical colour selec- 
tion criteria to exclude as many of the background galax- 
ies as possible, and then correct for background con- 
tamination among those galaxies with colours consistent 
with the spectrosc o pically -confirmed supercluster popula- 
tion. ICortese et al.l (|2008al ) found that Coma cluster mem- 
bers could be efficiently distinguished from most back- 
ground galaxies in the g — i versus NUV — i observed 
colour-colour plot. In the left panel of Fig. \TU\ we show 
the B — R versus NUV — R observed colour-colour plot 
used to produce an analogous separation between super- 
cluster members and the background population, while in 
the right panel we show the B — R/R colour- magnitude re- 
lation. The spectroscopically-confirmed supercluster mem- 
bers (large green symbols) lie within a well-defined region of 
the left-hand plot, while spectroscopically-confirmed back- 
ground galaxies (vel>17 000km s _1 ; large red symbols) gen- 
erally have redder B — R colours than supercluster mem- 
bers with the same NUV — R colour, although there re- 
mains significant overlap. We thus identify by eye suitable 
selection criteria to produce a compromise between simulta- 
neously minimize the number of supercluster members ex- 
cluded and the number of background galaxies included as 



indicated by the red lines. We also apply a second criterion, 
excluding galaxies lie above the cluster red sequence (right 
panel), where among the galaxies with redshift information, 
all are found to be background galaxies. Using this method 
we are able to exclude the bulk of the remaining galaxies 
for which we have no redshift information (blue points), al- 
though as apparent from the number of spectroscopically- 
confirmed background galaxies satisfying our selection cri- 
teria, there remain a significant level of background contam- 
ination among our colour-selected supercluster population. 

Figure [TOl shows that even with these colour cuts, a sig- 
nificant fraction of the remaining galaxies are background 
contaminants. We attempt to account for these statisti- 
cally, by measuring the real fraction, fsc(R), of superclus- 
ter members among those galaxies with redshifts satisfy- 
ing our UV-optical colour selection criteria for each bin of 
unit i?-band magnitude. This fraction declines slowly from 
1 for i?<15 to 0.80±0.04 for 17<ft<18 and 0.65±0.10 for 
18<7?<19. Following the completeness-correction method of 
de Propris et al. (2003), we assume that within each _R-band 
magnitude bin the fraction of galaxies that are supercluster 
members is the same in the spectroscopic subsample as in 
the photometric one. This should be true if the spectroscopic 
targets were selected according only to their i?-band mag- 
nitude and not their colour, which is largely true for each 
of the redshift surveys that contribute to our spectroscopic 
sample. It is important to note that we cannot apply a sim- 
ple completeness correction to galaxies as a function of their 
UV magnitude, as the spectroscopic sample within a given 
UV magnitude bin is heavily biased towards the quiescent 
population with NUV — R > 5. 

The FUV (blue) and NUV (red) galaxy luminosity func- 
tions obtained using this second method is shown by the 
solid symbols in Fig. 1111 The requirement for optical data to 
perform the UV-optical colour cuts means that these LFs are 
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Figure 11. The FUV (blue) and NUV (red) LFs for the whole 
2.8334 deg 2 region covered by both GALEX and WFCAM K- 
band imaging obtained using: (i) the statistical subtraction of field 
galaxy method (open circles); and (ii) UV-optical colour selection 
(solid circles). The dashed and solid lines indicate the best- fitting 
Schechter functions to the data using the two methods. 



based on the 2.26 deg 2 region covered by both the GALEX 
and SOS optical data, rather than the slightly larger region 
covered by WFCAM K-haxiA imaging (see Fig.[T]). The error 
bars include Poisson uncertainties as well as uncertainties 
in fsc(R)- We fit single Schechter functions (solid curves) 
via a x 2 - mmmi i za tion routine, obtaining best-fit values 
M^ uy =-18.27to\tl and aftfv=-1.505±0.12 for the FUV 
LF, and Mn UV =-18.76±0.31 and « NU v=-L505±0.085 for 
the NUV LF. Again the single Schechter fit provides a good 
description for the NUV and FUV LFs down to the survey 
limit. The best-fit parameters and their errors, calculated 
as the range of solutions within 1.0 of the minimum \ 2 are 
presented in Table 1 an d compared to previous determina- 
tions for the local field llWvder et al,ll2005l: iBudavari et"aD 
120051 ) , the Coma cluster~([Cortese et alj l2008al ) and Abell 
1367 (|Cortese et alj 120051 ). The errors in a and M* are 
highly correlated, and Fig. ll2l shows the joint 1, 2 and 3<r er- 
ror contours projected into the M* —a plane (solid curves) , 
while the dashed contours indicate X 2 ~Xmin~ LO. 

The ultraviolet LFs produced now are ~2o shallower 
than those obtained by the first method. We note that these 
two methods for the background correction are subject to 
different systematic errors, which are not fully quantified in 
the error bars. The difference between these results is likely a 
reflection of these systematics, and suggests that either the 
statistical subtraction of field galaxies is underestimating 
the background galaxies or the colour selection is missing a 
fraction of real members. We note that a number of previous 
claims of steep faint-end slopes in clusters obtained via the 
statistical subtraction method have been shown to be incon- 
sistent with shallo wer slopes obtained via deep spectroscopic 
follow-up surveys (R mcs fc Gelledl200sl ; IPennv fc Conselicd 




Figure 12. The confidence intervals of the Schechter function 
fits projected into the M*—a plane for the FUV (blue) and NUV 
(red) LFs shown in Fig. 1111 The solid contours delineate x 2 ~ 
Xmin = 2.30, 6.17 and 11.8, corresponding to the joint 1, 2 and 
3(7 uncertainties on M* and a, while the dashed contours indicate 
X 2 — Xmin = 1-0. The cor r espon ding local field galaxy FUV and 
NUV LFs of lWvder etHl j2005ll are also shown as blue and red 
points respectively with error bars. 



Table 2. Best-fitting parameters for the NUV and FUV LFs 



Band 


Sample 


Schechter parameters 
M* a 


NUV 
NUV 
NUV 
NUV 
NUV 


Shapley 
Coma^ 
Abell 1367^, 
Local field c 
0.07<z<0.13 field^ 


-18.76±0.31 
-18.50±0.50 
-19.77±0.50 
-18.23±0.11 
-18.54±0.15 


-1.505±0.085 

-1 77 +0.16 
L -' '-0.13 

-1.64±0.21 

-1.16±0.07 
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2008, e.g.). In these cases the steep faint-end slopes appear 
due largely to the contribution of background galaxies red- 
der than the C-M relation, and which are robustly removed 
by our UV-optical colours selection method. The statistical 
subtraction method is particularly suspectible to systematic 
errors due to photometric calibration offsets or corrections 
for Galactic extinction (which is large in our case). As field 
galaxy counts increase much more steeply than cluster mem- 
ber counts, small systematic uncertainties in background 
subtraction can produce large uncertainties in the numbers 
of faint cluster galaxies. For both methods used, the NUV 
and FUV LFs of the Shapley supercluster gal axies are con- 
sistent within errors with those produced by ICortese et alj 
(|2008a|) for the Coma cluster, while the My V values for Abell 
1367 appear a magnitude brighter, due to the presence in 
A 1367 of a few galaxies with enhanced star formation, in- 
cluding an ultra-luminous UV galaxy. These supercluster 
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Figure 13. The FUV LFs for the passive (NUV-R>4.5; red 
symbols) and star-forming (NUV — i?^4.5; blue symbols) super- 
cluster galaxy sub-populations. The solid lines indicate the best- 
fitting Schechter functions to the data. 



ultraviole t LFs are all significantly steeper than those ob- 
tained bv lWvder et al.l ([2005) for the local field galaxy pop- 
ulation (qjv(/v=-1-16±0.07, «Ft/v = -1.22±0.07; red/blue 
points with error bars in Fig. 11 2 p . 

Figure [13] examines the relative contributions of 
passive (NUV— J?>4.5; red symbols) and star-forming 
(NUV— R^4.5; blue symbols) galaxies to the far-UV lu- 
minosity function. We obtain best-fit Schechter func- 
tions for the star-forming supercluster galaxy popu- 



lations of Mftjv=— 18.12 



+0.37 



otFW =— 1.38_ 



and 



1.36^n'j ii consistent ( within 



M^ vv =-18.76t° ii, ajvuv=-1.36t°;i> 
the la errors) with those obtained by IWvder et al.l (1 20051 ) 
and iBudavari etall l|2005h for local field galaxies. The con- 
tribution of the passive galaxies to the far-ultraviolet LF 



is negligible at the bright end (Mfl 



r<- 



17), but becomes 



increasingly important at fainter magnitudes. 

It is important to note that the FUV emission from 
these "passive" galaxies is unlikely to be due to low-levels 
of recent or ongoing star-formation . In a HST/ ACS far- 
UV imaging survey ISalim fc Richl (|2010T ) also find that 
NUV— r~4.5 separates the early-type galaxy population 
into passive red sequence galaxies and those with evidence 
of extended residual star-formation in the form of rings 
or spiral arms. The FUV emission of these galaxies ap- 
pears extended in the GALEX image s, qualitative l y rul - 
ing out a dominant AGN contribution. lHaines et al.l (2008) 
show that there is very little contamination of the NUV —r 
red sequence by galaxies spectroscopically classified as star- 
forming (with Ha equivalent widths >2A) in the local uni- 
verse, while we find just one galaxy classified as star- forming 
among the ~100 NUV—R red seq uence galaxie s withi n the 
AAOmega spectroscopic sample of ISmith et al l (|2007l ). Us- 
ing the Rose Can index. [Smith et al. ( 20091 ) also excluded 
frosting by small mass fractions (1-5%) of recent star- 
formation (0.5-1 Gyr old) for the vast majority of their qui- 
escent galaxies in the AAOmega sample, which make up our 



NUV—R red sequence. Instead the near-UV flux in these 
quiescent galaxies is d ominated by hot main-sequen c e star s 
close to the turnoff (|Dorman. O'Connell fc Rood! 120031 ). 
while the far-UV radiation is thought to be produced mainly 
by old (>10Gyr) low-mass, helium burning stars in ex- 
treme horizonal- branch and subsequent phases of evolution 
jO'Connellllr999l ). 



4 MID/FAR-INFRARED LUMINOSITY 
FUNCTIONS 

We calculate the 24/im luminosity function for the 2.25 deg 2 
region covered by 24/im, optical and GALEX UV data. It is 
determined using the same completeness-correction method 
described in § 13.21 identifying possible supercluster mem- 
bers from their UV-optical colours, and statistically ac- 
counting for background galaxies using the existing spec- 
troscopic information. To account for the incompleteness of 
the 24/im data, we use the inverse of the completeness func- 
tion (shown in Figure [3]) as weighting factors when we cal- 
culate the number counts in each bin in 24/tm flux. The re- 
sulting 24/mi luminosity function, shown by the solid green 
points in Fig. 1141 covers a factor 300 in luminosity, extend- 
ing down to our survey limit of 0.35 mjy. The open sym- 
bols indicate the contribution just from the 360 galaxies 
with /24>0.35mJy and redshifts placing them in the Shap- 
ley supercluster, which represents >80% of galaxies with 
/24>8mJy, falling to ~40% at the completeness limit. We 
have redshifts for all /24>26mJy galaxies, and so are not 
missing any LIRGs in the supercluster, irrespective of their 
UV-optical colours. 

To estimate the total infrared luminosity, Lir(8- 
1000/im ), from ou r 24mn fluxes, we adopt the calibration 
used by IBai et alj (|2009l ) for the cluster A 3266 at z=0.06 
(after accounting for the small redshift difference between 
the Shapley supercluster and A 3266), which was derived 
by compariso n to the samp l e of st ar forming galaxy SEDs 
developed by iDale fc Heloul (|2002l ). These SEDs are based 
on IRAS and Infrared Space Observatory obser vations of 69 
normal star- forming galaxies (Dale et al. 2 00 ll ) and further 
calibrated with submillimeter observations. The SEDs are 
luminosity dependent, such that the 24/im contribution be- 
comes increasingly important with Lir, resulting in the lin- 
ear relation log L_f_R(ergs _1 )=42.91 + 0.86x log/24(mJy) at 
the distance of Shapley (shown along top axis of Fig. I14[) . 

We obtain a best-fit Schechter function of 
/ 2 V m =24.83t^mJy and a(24 M m)=-1.425t °^, 
shown by the dark green curve in Fig. 1141 Adopting 
the Lm—f2i relation derived above, these correspond to 
log(£W£©)=10.52±°;o1 and a(/i?)=-1.49±0.04. These 
values ar e fully consistent with both those obtained by 
IBai et all (|2006l ) for Co ma (log(£r R /Lp)) =10.48ln1? and 
a(7\R)=-1.49±0.11) and IBai et alj (|2009l ) for Abell 3266 
(log(LJ H /L Q )=10.49lg|? for a fixed a = -1.41). 

We now compare our supercluster 24/im galaxy LF with 
that for local field galaxies by overlaying a s blue squares 
the 24 /im luminosity function obtained by iMarleau et al.l 
|2007l ) for 0<2<0.25 star-forming galaxies in the Spitzer 
First Look Survey. We find the shape of the 24/xm LF of 
Shapley supercluster galaxies to be fully consistent with 
that obtained for local field galaxies, and support the as- 
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Figure 14. 24/im luminosity function for the Shapley superclus- 
ter (solid green symbols). The contribution due to spectroscopi- 
cally confirmed supercluster members is indicated by open sym- 
bols. Blue squares represent the field 24pm LF of iMarleau et au 
(2007) for < z < 0.25 star-forming galaxies in the Spitzer First 
Look Survey. Along the top axis we show the corresponding to- 
tal infrared luminosity (L/fl[8-1000/im]), which for ease of com- 
parison w i th the prev ious 24/xm cluster luminosity functions of 
iBai et al.l ll2006l. l2009h we deri ve directly from the L IR -f 2 4 re- 
lation used bv lBai et alj 11200911 for the cluster A 3266 at 2=0.06, 
after correcting for the change in expected 24pm flux levels be- 
tween Shapley and A 3266. 



sertion of IBai et al] (|200Sh that there is no environmental 
dependence on the sh ape of the 24pm luminosity function. 
IMarleau et"ail l|2007t) find evidence for an upturn in the 
galaxy counts below iAL 24Mm ~10 8 L© (Ljjj~5x10 42 ergs" 1 ), 
particularly once they include the contribution of galax- 
ies with only photometric redshift estimates, and we find 
marginal evidence for this upturn, albeit only in the faintest 
1-2 luminosity bins, where our spectroscopic complete- 
ness is lowes t and photometric uncertainties are greatest. 
IWestra et all (|2010h also find a similar upturn in the local 
(0.01<z<0.10) extinction-corrected Ha luminosity function, 
below L(//a)~10 40 ergs- 1 (SFR~O.1M yr" 1 ). 

We do not present an estimate of the 24pm supercluster 
LF using the background subtraction method, even though 
robust measure s of 24pm g a laxy counts down to our sur- 
vey limit exist. IShupe et al. (|200&f ) presented 24pm galaxy 
counts from the six SWIRE fields, covering in total 64deg 2 , 
but showed that even with these large fields (5-11 deg 2 ) con- 
taining tens of thousands of galaxies, the field-to-field dif- 
ferences were of the order 10%, much larger than could be 
accountable to Poisson uncertai nties. In the sub-mjy range 
(0.3-1 mJv UShupe et al.l (|2008l ) find that only 10-20% of 
24pm extragalactic sources are at «<0.3, and so the field 
24pm count variations, significant even on these 3 degree 
scales, may be due to large-scale structures at z>0.5. For 
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Figure 15. 70pm luminosity function for the Shapley super- 
cluster (solid green symbols). The contribution due to spectro- 
scopically confirmed supercluster members is indicated by open 
symbols. The best fitting Schechter function is indicated by the 
solid green curve and the associated parameters and uncertain- 
ties indicated. The blue squares indicate the local field luminosit y 
function of IRAS galaxies fromfyVang & Rowan- Robi nson] ||2010|) , 
while the blue dashed curve indicat es the analytic f o rm of the lo- 
cal field IR luminosity function of iTakeuchi et al I ll2003r> , both 
obtained from the PSCz survey. 



our faintest 24pm luminosity bin, the field-to-field variance 
seen in the SWIRE galaxy counts (D. Shupe, private com- 
munication) corresponds to an uncertainty of ^400, double 
the estimated cluster galaxy counts from Fig. [14] The huge 
uncertainties due to cosmic variance can be understood as a 
consequence of just ^5% of the 24pm sub-mjy extragalac- 
tic sources in the Shapley field actually belonging to the 
supercluster. Note that the uncertainies obtained using the 
completeness-correction method as indicated by the error 
bars in Fig. [14] are much smaller, due to a combination of 
the extensive spectroscopic information, and the ability of 
the UV-optical colour selection to exclude the bulk of the 
background 24pm sources. 

In Fig. [TS] we present the corresponding 70pm lumi- 
nosity function of Shapley supercluster galaxies (solid green 
symbols) using again the completeness-correction method of 
§ 13.21 As before, the open symbols indicate the contribution 
spectroscopically-confirmed supercluster members. Along 
the top axis we show the correspon ding bolome t ric inf rared 
luminosites based on the SEDs of iRieke et al.l (|2009l ). We 
fit the 70pm luminosity function by a single Schechter func- 
tion as shown by the green curve, obtaining best-fit parame- 
ters of^^^^lQSt^fmJy and a(70pm)=-1.40±o|5. Using 
the lRieke et al" (2009) SEDs, we convert this to a bolomet- 
ric infrared luminosity, obtaining log(L| fl /Z/Q) = 10.44lQ 22 . 
This is the first measurement of the 70pm luminosity func- 
tion of cluster galaxies with Spitzer, and appears fully con- 
sistent with that obtained at 24pm, albeit with large un- 
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certainties in the Schechter parameters due to the rela- 
tively shallow depth. We compare the supercluster 70/im 
with the 60/xm local field galaxy lumino sity function ob- 
tained by [Wang fe Rowan-Robinson] (|2010h from the all-sky 
IRAS Point Source Redshift Survey (blue squares), assum- 
ing a MIPS-70 / jm/IR AS-60/jm flux ratio of 1.56 based on 
the lDale et"afl l|200ll ) infrared SED with a=2.5, redshifted 
to z = 0.048. This model SED best describes the /70//24 
mid-infrared colours of the majority of our 70^tm-bright su- 
percluster population (Paper III). The blue dashed line in- 
dicates the best-fit analytic f unctio n for the local 60/xm LF 
obtained by iTakeuchi et all {2003). Again, we find no sig- 
nificant difference between the 70^im luminosity function of 
Shapley supercluster galaxies and that obtained for local 
field galaxies. 



5 COMPARISON OF THE UV AND MIR 
EMISSION FROM CLUSTER GALAXIES 

The infrared-to-ultraviolet ratio is a rough measure of the 
amount of extinction at ultraviolet wavelengths of young 
stellar populations (<10 8 yr), with the caveat that the ob- 
served ultraviolet and infrared emission do not come from 
exactly the same regions within a star-forming region or 
galaxy, or be produced by recent ly formed stars of the same 
ages/masses (Calzetti et alj2005f ). The main problem affect- 
ing UV and optical-based SFR indicators is dust obscura- 
tion, which can require corrections of the order 1 0-100 in th e 
SFR estimator, particularly in the ultraviolet (Bell 2002). 
This correction is also sensitive to the metal content, star- 
formation history and the relative distribution and geometry 
of the interstellar dust with respect to the s t ars producing 
the ul t raviolet emission ([Calzetti et al.l l 2005l ; ICortese et al.l 
l2008al ; iBoauien et alj|2009h . 

In Figure \W\ we show the infrared-to-ultraviolet ra- 
tio as a function of Tf-band magnitude for those galaxies 
spectroscopically-confirmed as Shapley supercluster mem- 
bers, estimating the total infrared luminosity (Ltir) for 
each gala xy from thei r 24mn fluxes using the infrared SED 
models of lRieke et all (|2009). The right-hand axis shows the 
resultant level of attenuation in the FUV estimated from 
the infrared-to-ultraviol et ratio ba s ed on the prescription 
described in Eq. 2 from lBuat et al.l l|2005l ). This is valid for 
star-formation histories in which the bulk of dust emission 
is related to the ultraviolet absorption, as for most nor- 
mal star-forming galaxies, and any range of configuration 
of dust /star geometry, metallicity or atte nuation law shoul d 
affect the relation by less than 20 per cent (|Buat et al]|2005t ). 
In cases where a signifcant fraction of the dust heating comes 
from evolved stars, such as early-type galaxies, the level of 
attenuation in the FUV should be significantly lower fo r 
the same infrared-to-ultraviolet ratio ICortese et al.l (j2P08b). 
The colours of each symbol indicate the /24//K flux ratio, 
which can be considered a proxy for the specific-SFR of the 
galaxy, from orange/yellow (0.15-1), through green (1-8; ly- 
ing along the star-forming sequence, consistent with being 
normal spiral galaxies) to blue (>8; above the main star- 
forming sequence). In paper III, we give a full discussion of 
how the observed /24//K flux ratio relates to the specific- 
SFR and the bimodality between star-forming and passive 
galaxies. In the case of passively-evolving galaxies neither 
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Figure 16. Infrared-to-ultraviolet ratio as a function of X-band 
magnitude for star-forming galaxies (/24//K >0.15) in the Shap- 
ley supercluster. Symbol are coloured according to their /24//K 
ratio from yellow/orange (<0.8) to blue (>8). Lower limits are 
shown for galaxies fainter than our nominal far-ultraviolet com- 
pleteness limit of mpuv='22.5. The diagonal line indicates a trend 
of the form Ltir/ LpuyocLx ■ The far-ultraviolet ext inction pre- 
script ion used for the right-hand axis comes from iBuat et al.l 
(|2005h . 



the FUV nor the 24/im emission is believed to be primarily 
produced by star-formation, but rather evolve d stars (see 
Paper III; IO'ConneIilll999l ; iBressan et alj|2006j ). and so we 
show only galaxies identified as star-forming according to 
having /24//K >0.15. 

The infrared-to-ultraviolet ratio correlates strongly 
with the if -band luminosity, increasing from ~1 at K^16 
(.M~10 9 Mq) to ~100 for the most massive galaxies, consis- 
tent with an overall trend of the form Ltir/ LfuvocLk (di- 
agonal line) . The trend for more massive galaxies to be more 
heavily attenuated in the ultraviolet than low-ma s s sys - 
tem s was previously d escribed by both lCortese et all (|2006t ) 
and lDale et al.l (|2007l l. although both these authors obtain 
somewhat shallower trends than found here. However, at 
fixed A"-band luminosity, there remains a significant scat- 
ter, which also increases with A"-band luminosity. These 
tre nds are consistent also w ith t hose obtained f or A (Hq) 
by iBrinchmann et all (|2004T ) and iGarn fc Best] (|2010h for 
starforming galaxies in the SDSS survey, as estimated from 
the Balmer decrement (Ha/H/3). The primary driver be- 
hind the increase in extinction with stellar mass is thought 
to be a combination of the greater dust content of more 
massive galaxies (as would be produced assuming a con- 
stant dust gas fraction for example) and the mass-metallicity 
relation, such that more massive (a nd hence metal-rich) 
galax ies have higher dust-to-gas ratios (iMunoz-Mateos et al.l 
l2009t ). Finally, at fixed 7f-band luminosity, the infrared-to- 
ultraviolet ratio increases with /24//K , i.e. specific-SFR, the 
blue symbols appearing located systematically above the or- 
ange/yellow ones. 
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Figure 17. Comparison of the FUV and 24/im luminosity func- 
tions in terms of the derived SFRs, assuming the UV a nd MIR 
components of the SFR calibration of lLerov et al.l (2008). 

5.1 The global contributions of obscured and 
unobscured star-formation 

As we have just seen, the star formation in many of the 
supercluster galaxies is heavily obscured, with more of the 
energy from young stars output in the infrared than the 
ultraviolet. How does this average out over the entire super- 
cluster popula tion? Using the ul traviolet and 24pm-based 
calibrations of lLerov et alj (|2008l ) 

SFR(M yr _1 ) = 0.68x 10~ 28 L„(FUV) [erg s^Hz" 1 ] 

+ 2.141^9 xlO- 43 L(24 M m)[ergs- 1 ], (1) 

(which assumes a Kroupa IMF), we can directly compare 
the observed FUV and 24pm luminosity functions in terms 
of star-formation rates, the result of which is shown in Fig- 
ure !17l We can immediately see that for the galaxies with the 
highest SFRs, the estimates obtained via the mid-infrared 
are ~10x higher than those obtained from the uncorrected 
far-ultraviolet, and that for SFRs>l M yr -1 , the far-UV 
component can be neglected and the 24/xm LF can be used 
as a good estimator of the SFR distribution. At lower SFRs 
however (^lMgyr -1 ), the FUV component cannot be ne- 
glected, as it becomes equally important as the infrared com- 
ponent, comprising anywhere between 10-90 per cent of the 
emission from star-formation in these galaxies. This trend 
for the dust attenuat ion to increase wi th SFR is consistent 
with that observed bv lBuat et all (|2007l) for local field galax- 
ies, and suggests that for surveys limited to star formation 
rates >1 Mq yr -1 , mid-infrared observations alone should be 
sufficient, encompassing >90 per cent of the emission from 
ongoing star-formation. 

Over the entire SOS region, we estimate a total super- 
cluster SFR of 327teo 2 M y r ~\ of which 264tJ° 2 M yr" 1 
(^80 per cent) is obscured (based on the 24pm calibra- 



tion) and just 63±3MQyr _1 (~20 per cent) is unobscured 
and emitted in the form of ultraviolet continuum. For the 
mid-infrared contribution we excluded those galaxies with 
/24//K <0.15 for which we assume the 24pm emission comes 
from evolved stars give n that all such ga laxies within the 
spectroscopic sample of[Smith et all (2007) are classified as 
passive based on their lack of Ha emission (Paper III), and 
have colours similar to early-type galaxies whose infrared 
emission has been shown to come from the dusty circum- 
stellar envelopes of mass- losing AGB stars (|Bressan et al.l 
120061 ; IClemens et alj|2009h . 

The infrared error is entirely dominated b y the 30 per 
cent uncertainty in the SFR calibration of iLerov et al.l 
(2008), which is comparable to the differences among the 
most widely used 24pm SFR calibrations in the liter- 
ature (for a comparison of these see e.g. ICalzetti et al. 
l2010h . We note that the 24pm calibration o f ILerov et all 
l|2008h is consistent to that of iRieke et all (|2009h which 
is SFR( M a yr- 1 )=2.04x lQ- 43 L(24Atm) [ ergs~ 1 1, as well as 
those of lWu et al.l (|2005h and lZhu et all (|2008T ). all of which 
are base d on global galaxy measurements. Using instead the 
SEDs of lRieke et al.l l|2009l ) to estimate the total infrared lu- 
minosi ty from the 24jxm fluxes and the infrared SFR calibra- 
tion of lBuat et al.1 (1200 8) . we obtain SFRs that are systemat- 
ically 10-20 per cent lower, once we account for the expected 
30 per cent contribution to the infrared flux from evolved 
stars (|Bellll2003l ; llglesias-Paramo et al.ll2004h . although we 
note that this relative contribution is likely to depend on the 
part icular star-formation histories of the gala xies involved 
(e.g. ICortese et al.|[2008bl ; ICalzetti et al.ll2010h . Even better 
consistenc y is found if we use instead the Ltir-SFK cali- 
bration of lBelll (|2003l ). producing SFRs within 10 per cent of 
those based on Eq. [T]for ~10 10 — 10 11 L© after accounting 
for the different IMF used. 

The remaining sources of uncertainty contribute to an 
uncertainty in the SFR/_r of just 11 Mq yr -1 . These include 
the 4% absolute calibrat ion uncertainty of MIPS at 24pm 
l|Engelbracht et all 120071 ). the statistical uncertainty aris- 
ing from including the photometrically-selected superclus- 
ter population, and the effects of varying the cut-off in the 
range 0.1</24//k<0.2. For the far-ultraviolet contribution 
we excluded those galaxies with NUV — R>4.5 (the emission 
again due to evolved stars; see § 3.2), and the error estimate 
includes the 0.05 mag uncertain ty in the FUV absolute cal- 
ibration (|Morrissev et al.l 120071 ) and statistical uncertainty 
from the photometrically-selected supercluster members. 

We believe the bulk of the mid-infrared emission (and 
hence our SFR budget) to be due to star formation, rather 
than AGN, whose contribution we estimate as of the order 
5-10 per cent. This we base on the observed tight FIR- 
radio correlation seen for the far-IR bright cluster members 
(Paper III) when comparing both the 2 4pm and 70pm fluxes 
with the published 1.4 GHz catalogue of lMillerl (120051 ). More- 
over, these galaxies have FIR-radio colours more consis- 
tent with dust heated by low-intensity star formation rather 
than AGN, and have extended 24pm emission. However in 
Paper III, we find four 24pm-bright sources whose emis- 
sion remains unresolved by MIPS, and whose FIR-radio 
colours are more consistent with dust heated b y AGN. Two 
of the se are within the spectroscopic sample of ISmith et al.l 
l|2007h and are classified as AGN, and between then con- 
tribute ^15Moyr -1 to the cluster SFR budget. Unfortu- 
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nately, we do have Spitzer/IRAG photometry, to identify 
IR-bri ght AGN from the ir power-law spectra in the 3-10/im 
range jstern et alJl2005l ). making a complete census of the 
AGN population impossible. There is however partial X-ray 
coverage of the SCC via a mosaic of 10 XMM images, each 
with exposure times of 20-40 ksec, covering A 3562, A 3558, 
SC 1329-313 and SC 1327-312 and most of the connecting 
filamentary strucutre, allowing us to identify AGN as X-ray 
point sources. Based on a sample of 202 spectroscopic SSC 
members detected at 24/im, comprising 52 per cent of the 
24/im flux associated to SSC galaxies, lying within the XMM 
mosaic, we estimate that IO.6I42 per cent of the 24/tm flux 
comes from X-ray loud AGN. We note that in infrared-bright 
sources, often both AGN and starburst activity are ongoing 
concurrently, and that the contribution from the AGN com- 
ponent (typically 5-50 p er cent) can only be resolved v ia in- 
frared spectroscopy (e.g. Goulding fc Alexander! [2009). and 
so we believe that this estimate of ~10per cent AGN con- 
tamination represents an upper limit. 



6 DISCUSSION 

The luminosity function is one of the most basic and funda- 
mental properties of the galaxy population, and its form 
is shaped by many of the processes of galaxy formation 
(| Benson et al.|[2003l ). The optical and near- infrared luminos- 
ity fun ctions of fi e ld an d cluster galaxies are well described 
by the ISchechterl (| 19761 ) function, which has a sound theo- 
retical basis as the expected mass function of galaxies built 
up by gravitational hierarchical assembly, although it seems 
that additional feedback processes from AGN and super- 
novae are required to fit the observed sharp cut off at bright 
lumin ositi es and the shallow fain t -end slope (1 Bower et al.l 
120061 ). As lYun. Reddv fc Condor] l|200ll) argue, if star for- 
mation activity is an integral part of galaxy evolution and 
reflects the nature of the host galaxies, then the Schechter 
function should also offer as good a description for the ultra- 
violet, far-infrared and radio LFs as it does for the optical 
and near-infrared LFs. 

In Sections [3] and [4] we have obtained the ultravi- 
olet and far-infrared galaxy luminosity functions for the 
Shapley supercluster, complementing the existing optical 
jMercurio et al.l 120061 ) and near-infrared l|Merluzzi et al.l 
1201(f ) LFs for the same region. In the ultraviolet, using 
new GALEX imaging, we find that the NUV (A e =2310A) 
and FUV (A e =1510A) LFs are both well described by sin- 
gle Schechter functions down to Mw=— 14 (~M*+4.5) 
with Mfe/v=-18.27±0.41, aFi/v=-1.505±0.12 for the far- 
UV, and Mjvuv=-18.76±0.31, ajvc/v=-l-505±0.085 for 
the near-ultraviolet. These LFs are essentially identical once 
the typical ultraviolet colours of star-forming galaxies found 
in the supercluster having FUV— NUV^0.5 (FxcxX^ 1 ) are 
taken into account. These ultraviolet LFs of the Shapley su- 
percluster core are fully consiste nt with those obtained by 
ICortese et all (|2003l . l2005l . r2008al ) for Virgo, Coma and Abell 
1367. 

The faint-end slopes of local cluster ultraviolet LFs 
(q~— 1.5) are all observed to be significantly steeper (at 
the ^2a lev e l) th an th e local field ga l axy U V LFs of 
IWvder et af] (|2005l 'l and iBudavari etall (|2005l ), implying 
that the shape of the LF depends on environment (clus- 



ter versus field). This appears at least partly due to the 
morphology-density or star formation-density relations, as 
Fig. \n\ shows that the blue, star-forming (NUV-R<4.5) 
and quiescent (NUV— i?>4.5) galaxies contribute to the 
LFs in quite different ways. The blue, star-forming galax- 
ies dominate at the bright-end, and have a faint-end slope 
aF[7v = 1.38io '135 consistent within errors with the local 
field galaxy FUV LF (albeit still marginally steeper). The 
contribution of the passive galaxies to the far-ultraviolet 
LF is negligible at the bright end (Mfuv<— 17), but be- 
comes increasingly important at fainter magnitudes. This 
pile-up of quiescent galaxies in the last three magnitude 
bins (Mfuv> — 15) appears the primary cause of the dis- 
crepancy in the faint-end slopes between cluster and field 
galaxie s. This effect is seen a lso in the Coma cluster ultravio- 
let L F JCortese et al|2008al ) , and also at optical wavelengths 
(e.g. iMercurio et al.ll2006l ). A second contributing factor to 
the steeper observered cluster UV LFs could be simply the 
steeper stellar mass fu nctions (or if -band L Fs) of cluster 
galaxies, as observed bv lMerluzzi et all (|2010l) for this same 
Shapley supercluster core region, finding olk = — 1.42±0.03 
as o pposed to cxk=— 1.16i:0.04 for the local field if-band 
LF \ Jones et al.ll2006l ). Given that the FUV emission from 
the "passive" galaxies is unlikely to be due to star-form ation 
jO'Connelll 1 19991 : iDorman, O'Connell fc Roodl 120031 ). the 
relative consistency between the FUV LFs of cluster star- 
forming galaxies and the field, suggests little environmental 
dependance in the unobscured SFRs of star-forming galax- 
ies. 

In §2] using new panoramic Spitzer/MlPS 24/im imag- 
ing, we determined the infrared galaxy LF for the Shap- 
ley supercluster core (Fig. I14|) . finding it well described by 
a single Schechter function with log(Z/| fl /Z/Q) = 10.52+Q'[]g 
and a(IR)—— 1.49±0.04. Thes e values a r e fully consistent 
with both those obtained by iBai et all |2006) for Coma 
(log(L* f R ALb ) = 10.481^1? and a(7J?)=-1.49±0.11) and 
IBai et al.l l|2009l ) for Abell 3266 (log(Ll fl /L )=lO.49^i? 
for a fixed a = —1.41). Although we do not extend to quite 
as low luminosities as IBai et all (|2009l ) our larger sample 
size allows us to derive stronger constraints on the 24/im 
luminosity function of local cluster galaxies. 

The Spitzer/MIPS 24/im band is some way from the 
typical peak wavelength (~100/im) of the far-infrared emis- 
sion from galaxies, requiring uncertain extrapolations in es- 
timating Lib. f rom the 24/xm flux su ch as via the m o del in - 
frared SEDs of iDale fc Heloul ||2002D or iRieke et ail l|2009h . 
The Spitzer/MIPS 70/im filter, being closer to the far- 
infrared peak affords more reliable Lib. estimates, and so we 
have redetermined the infrared LF, based on our available 
70/im photometry of the Shapley supercluster core (Fig. [To)) . 
A best-fit Schechter function with log(L| J j/L Q ) = 10.44l° ^ 
and a(77?)=-1.40j:^^ is obtained, fully consistent with 
the 24/tm LF, albeit with large uncertainties due to the 
relatively shallow depth of our 70/im imaging. This rep- 
resents the first measurement of the 70/im galaxy LF of 
a local cluster with Spitzer, and should prove a useful 
local benchmark to follow the evolution in the infrared 
properties of cluster ga l axies w it h ongo i ng and future sur- 



veys (e.g. Haines et al. 2009a 



iBraglia et all l2010l : iRawle et al l 



2010al : ISmith et alj 
20ld : IChung et all 



2010; 
B01C ) 



with Spitzer and particularly Herschel whose PACS instru- 
ment covers rest-frame 70/im beyond z~l. 
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iBai et all (120061 ) found that the 24/im galaxy LF of 
Coma was consistent with the infrared LFs of local field 
galaxies, based on comparison to the 12/x m LF derived from 
the all-sky IRAS Faint Source Catalogue (|Rush et al.lll993l ^ 
or the 60/mi LF o f local galaxies in the IRAS Point Source 
Redshift Survey l|Takeuchi et a.1.1 l2003h . We also find the 
shape of both the 24/im and 70/im LFs of Shapley super- 
cluster galaxies to be full y consistent wi t h tha t obtained 
for local field galaxies by iMarleau et al ] J2007I) at 24/im 
(blue squares in Fig. [Til and [Takeuchi et all i|2003t ) and 
IWang fc Rowan- Robinsonl |2010f ) at 70mn (dashed line and 
blue s quares respective ly in Fig. I15[) . confirming the asser- 
tion of IBai et all ((2009) that there is no environmental de- 
pendence on the shape of the far-infrared luminosity func- 
tion. 

A probable significant consequence of this apparent 
complete agreement between the local cluster and field 
galaxy infrared LFs, is that the bulk of star-forming galaxies 
that make up the observed cluster infrared LFs have been 
recently accreted from the field and haven't had their star 
formation activity significantly affected by the cluster en- 
vir onment yet, based o n a similar argument to that used 
by iBaloeh et alJ ((20041 ) regarding the lack of environmen- 
tal dependence seen for the distribution in EW(Ha). If the 
mechanism which quenches star-formation in infalling galax- 
ies acts rapidly, then these galaxies rapidly drop out of the 
sample, and no longer contribute to the far-infrared LF (ex- 
cept perhaps at the very lowest luminosities studied here). 
The primary change to the infrared LF is then a reduction 
in the overall normalization, but without affecting its shape 
(i.e. L* or a). If instead, the quenching of star- formation 
acts slowly (>1 Gyr), then a significant fraction of the clus- 
ter galaxies which constitute the infrared LF will be those in 
the process of being quenched, and have reduced (often by 
a factor 2 or more) SFRs and infrared luminosities in com- 
parison to the field population. The primary consequence 
of there existing a significant population of transforming 
galaxies would be that the faint-end slope steepens, while 
the Schechter function becomes an increasingly poor fit to 
the data. However, the shape of the infrared LF is not ex- 
pected to be a particularly sensitive indicator to the pres- 
ence of galaxies being quenched by en vironmental processes , 
based on similar arguments to those of lCortese et alJ (|2008aT ) 
about the ultraviolet luminosity function. 

The bulk of these infalling star-forming galaxies can- 
not be interlopers, who happen to have line-of-sight ve- 
locities within the redshift range of the cluster, as the 
normalization of the cluster infrared LF is still ~25-30x 
higher than the local field infrar ed LFs (|Marleau et al ] |2007l : 
IWang fc Rowan-Robinsonll2010h multiplied by the total vol- 
ume occupied by the cluster survey (~2 400 Mpc 3 ). We note 
that this higher normalization does not imply star formation 
is being triggered in the supercluster, but simply reflects the 
overall significant overdensity of galaxies in this volume. 

It is difficult to understand why if the field and cluster 
infrared LFs are consistent, the cluster FUV LF should be 
significantly steeper than that found in the field, as both re- 
flect the star-formation from the essentially same galaxies. 
One possibility is that the differences are only fully apparent 
at the very low SFRs traced by the GALEX data, and which 
are missed by the marginally (~2-3x) shallower 24/im data. 
Alternatively, the skewed FUV LF could reflect environmen- 



tal processes acting on just dwarf galaxies which are missed 
by our 24/xm survey, as their star formation is unobscured 
and emitted primarily in the FUV. 

One small caveat to the apparent agreement between 
cluster and field infrared LFs over the luminosity range 
10 9 — lO 11 !/© is the observed behaviour of the local field 
far-infrared LF, which for Lir^W 11 Lq can no longer 
be described by a single Schecter function. Instead of 
declining exponentially, the luminosity function is better 
described as a powe r- law (|Yun. Reddv fc Condon! l200ll : 
iTakeuchi et"afl 120031 ). or via an additional Schechter 
function to model the galaxy counts at Li r^W 11 Lq. 
Similarly, the 1.4 GHz radio galaxy luminosity function 
cannot be described by a single Schechter function, but 
requires two components, one representing the normal, 
late-type field galaxies and the other representing "mon- 
sters " powered by AG N dCpndon. Cotton Sz Broderickl 
120021 : iMauch fc Sadler! I2007T) which dom inat e at 
ii^GHz^lO^WHz" 1 . lYun. Reddv fc Condon! ((20011 ) 
thus model the far-infrared LF as a sum of two Schechter 
functions, one for normal field galaxies, and a second to 
model the high-luminosity excess composed of a star- 
burst /ULIRG population in which the SFRs (and hence 
FIR luminosities) are temporarily boosted by an order of 
magnitude. 

The apparent ability to fit the cluster far-IR LFs by 
a single Schechter function, could suggest that this second 
starburst component is missing from clusters. However, we 
note that the density of such luminous infrared galaxies in 
the field is <10" 5 Mpc- 3 , and given our survey volume of 
just ~2 400Mpc 3 we would expect to find <<1 such galaxy 
in the SSC. 



7 SUMMARY 

We have presented new panoramic Spitzer/MlPS mid- and 
far-infrared (MIR/FIR) and GALEX ultraviolet imaging of 
the most massive and dynamically active system in the lo- 
cal Universe, the Shapley supercluster at z=0.048, covering 
the five clusters which make up the supercluster core. Using 
existing spectroscopic data from 814 confirmed superclus- 
ter members, we have produced the largest complete cen- 
sus of star-formation (both obscured and unobscured) in lo- 
cal cluster galaxies to date, extending down to SFRs~0.02- 
0.05 Moyr -1 , providing a local benchmark for comparison to 
ongoing and future studies of cluster galaxies at higher red- 
shifts with Spitzer and Herschel. From the GALEX data we 
produced near-ultraviolet (NUV) and far-ultraviolet (FUV) 
luminosity functions, which were found to have steeper faint- 
end slopes (a— — 1.505) than the local field population at 
the ~2<r level, due largely to the contribution of massive, 
quiescent galaxies at Mfuv> — W- Using the Spitzer/MIPS 
24/im imaging, we determined the infrared galaxy LF of 
the Shapley supercluster core, finding it well described by a 
single Schechter function with log(L} R / Lq) — 10. 52^_gQg and 
a(IR)=— 1.49±0.04, fully consistent with those obtained for 
the Coma cluster and Abell 3266. We also presented the 
first 70/mi galaxy luminosity function of a local cluster with 
Spitzer, finding it to be consistent with that obtained at 
24/im. The shapes of the 24/im and 70/mi galaxy luminos- 
ity functions were also found to be indistinguishable from 



18 Haines et al. 



those of the local field population. This apparent lack of en- 
vironmental dependence for the shape of the FIR luminosity 
function suggests that the bulk of the star-forming galaxies 
that make up the observed cluster infrared LF have been re- 
cently accreted from the field and have yet to have their star 
formation activity significantly affected by the cluster envi- 
ronment. Combining the GALEX FUV and Spitzer 24/im 
photometry we estimate a global SFR of 327tgQ 2 M0yr _1 
over the whole supercluster core, of which just ~20 per cent 
is visible directly in the ultraviolet continuum and ~80 per 
cent of the energy from star formation is reprocessed by dust 
and emitted in the infrared. 
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